The Emergence of El-Nino as an Autonomous Component in the Climate Network 
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We construct and analyze a climate network which represents the interdependent structure of 
the climate in different geographical zones and find that the network responds in a unique way to 
El-Niho events. Analyzing the dynamics of the climate network shows that when El-Niho events 
begin, the El-Niho basin partially loses its influence on its surroundings. After typically three 
months, this influence is restored while the basin loses almost all dependence on its surroundings 
and becomes autonomous. The formation of an autonomous basin is the missing link to understand 
the seemingly contradicting phenomena of the afore-noticed weakening of the interdependencies in 
the climate network during El-Niho and the known impact of the anomalies inside the El-Nino basin 
on the global climate system. 
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It was recently suggested that climate fields such as 
temperature and geopotential height at a certain pres- 
sure level can be represented as a climate network 
In this network, different regions of the world are repre- 
sented as nodes which communicate by exchanging heat, 
material, and by direct forces. These interactions are 
represented by the links of the climate network. Interac- 
tions between two nodes may also exist due to processes 
which take place outside the atmospheric pressure level 
or through interactions with the ocean and the lands. 
Each link is quantified by a weight based on measures of 
similarity between the time series (e.g. correlations) of 
the corresponding individual nodes (see [2] for an experi- 
mental evidence for the relations between heat exchange 
and synchronized fluctuations of the temperature field). 

Recent studies show that many links in the cli- 

mate network break during El-Nino events. The cli- 
mate network contains several types of links, that have 
different levels of responsiveness to El-Nino. From the 
maps of Tsonis and Swanson [4], one can locate the re- 
sponsive nodes (the nodes attached to the most respon- 
sive links) to be in the pacific El-Nino Basin (ENB) 0|. 
These maps together with the tremendous impact of the 
El-Nino Southern Oscillations on world climate, suggest 
that ENB have a unique dynamical role in the dynamics 
of the climate network. Indeed ENB has unique topolog- 
ical properties. Its connectivity and clustering coefficient 
fields studied by Donges et. al. [7] can be distinguished 
from their surrounding by their particularly higher val- 
ues. Also, the betweenness centrality field [8] in ENB is 
very low. However, the dynamics related to the interac- 
tion of ENB with its surroundings and the origin of its 
unique features are still not known [9j. 

In this Letter we follow the dynamics of the climate 
network in time, between the years 1979 and 2009 where 
eight El-Nino events took place. We identify a cluster 



around ENB that shows a clear autonomous behavior 
during El-Nino epochs, and determine the dynamics of 
its interactions with the surroundings. We also find an 
epoch of a decreased influence of ENB on the surround- 
ings, typically three months before the emergence of the 
autonomous behavior. Our findings resolve the seem- 
ingly contradictory situation of decreased interactions of 
ENB with its surroundings on one hand, as explored in 
previous works [3|— [Hj] , and the known influence of ENB 
on world climate on the other hand. We find that only 
the influence of the surroundings on the ENB region is 
significantly reduced, creating, apparently, a dynamical 
autonomous source within this region which influences 
its surroundings. 

The adjacency matrix in our climate network analysis 
is based on a similarity measure between time series of 
temperatures (after removal of the annual trend) at the 
pressure level of lOOOmb covering the last 30 years. In 
the current work we pick 726 nodes from the Reanalysis 
II grid [lOj , such that the globe is covered approximately 
homogeneously. The similarity measure Wf r we use is de- 
fined as follows. First we identify the value of the highest 
peak of the absolute value of the cross covariance func- 
tion. Then we subtract the mean of the cross covariance 
function and divide the difference by its standard devi- 
ation (see [3|, l5| for further details). The indices I and 
r represent two nodes, and y is the beginning date of a 
snapshot of the network, measured over 365 days with 
extended period of 200 days of shifts needed for evalu- 
ating the cross covariance. The time shift of the highest 
peak of the cross covariance function from zero time shift 
is denoted 9f r . The sign of Of r stands for the dynamical 
ordering of I and r (and hence Q\ r = — v rl ). In directed 
graph theory terms, when d v lr > the link is regarded as 
outgoing from node I and incoming to node r. Until now, 
such ordering or direction were not considered for con- 
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structing climate networks. Our method, thus, enables 
to distinguish between in- and out- links. 

The adjacency matrix of a weighted directed climate 
network is defined : 



Af ir = (l-«J,, P )e(6>» r )w£ 

whereas Q(x) is the Heaviside function. 

The in- and out- weighted degrees are defined : 



(i) 



if 



of = E A l 



(2a) 



(2b) 



The / and O fields represent the level of the depen- 
dence of a node on its surrounding, and the level of its 
influence on the surrounding, respectively. 

The total weighted degree (also denoted "vertex 
strength" of a node I is defined as : 
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Fig. [Tal shows the time dependence of D for each node 
in the network. Two main observations can be clearly 
seen from this figure. First, weighted degrees of the nodes 
yield an extremely persistent quantity [12], and different 
geographical regions have typical values over time. A 
second notable pattern is the horizontal bright stripes 
that appear during El-Niho events. The second feature 
further supports the finding [345| that many links of the 
climate network break during El-Niho events. 

Inside the mid range of node indices, between 300 and 
500, there is a group of nodes, C, that have lower values 
at any time, and are specifically distinct during El-Niho 
(Fig. [Taj). These low values and their distinct response to 
El-Niho, as we will explicitly show, are mainly due to in- 
links. When studying the in- and out- links separately, 
we find that the spatial distribution of the in- weighted 
degrees If is much broader compared to the distribution 
of the out- weighted degrees Of, at any time y, typically 
by 10 to 20 percent. The 90-th percentile (statistics is 
collected over time) yields a deviation of 30 percents. As 
shown below, this difference between / and O is mainly 
contributed by the group C discussed above. In order to 
follow the different dynamical behavior of / and O, we 
apply the following approach. 

We consider only links that are related directly to the 
El-Niho dynamics by focusing on the group of nodes C 
that have extremely low weighted degrees during El-Niho 
events: 
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FIG. 1: Weighted degrees as a function of time (y axis) 
and the node index (x axis), (a) Total weighted degree 
Df. (b) The microscopic contributions to the weighted 
incoming degree of C, IjL. (c) The microscopic 
contributions to the weighted outgoing degrees, Oq. 
One should bare in mind that each point is compiled 

from records of 565 days : 365 days + 200 days of 
shifts. The representative point in all figures for each 
565 period is the beginning date of the period 



where T = 2150 is a threshold [13|, and the triangular 
brackets stand for averaging over all El-Niho events. The 
group C which includes 14 nodes is part of ENB, and is 
notably located in the eastern equatorial part of ENB, 
which is known to have a large scale upwelling of cold 
ocean water, yielding a "cold tongue" [6] which deforms 
during El-Niho (see Fig. [6a] for its location). 

We find that measuring the dynamics of the interac- 
tions of the C region with its surroundings yields a sensi- 
tive tool to quantify the responses of / and O to El-Niho 
events. One can define 1^ and O v c , the in- and out- 
weighted degree of C, respectively, in the same manner 
we defined the weighted degrees for nodes. We consider 
only links between pairs of nodes where one node be- 
longs to C and the other does not. Fig. [2] shows the time 
dependence of these two variables. As seen in the fig- 
ure, the responses of the in- and out- degrees to El-Niho 
seem to be anti-correlated. When the m-degree drops, 
the out-degree slightly increases. The interpretation is 
that the nodes inside C lose large part of their depen- 
dence on the surrounding, while slightly increase their 
influence outside, thus becoming significantly more au- 
tonomous during El-Niho. Note that the excess of out- 
links over in- links is evident also in normal times. The 
detailed space-dependent picture, however, is even richer, 
as discussed below. 

The spatial dependence of Iq and O y c is explored by 
what we call the "microscopic" contributions from each 
one of the 712 nodes that do not belong to C. These 
microscopic contributions can be visualized (Figs. fllEtcj) 
as they evolve in time. 
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FIG. 2: Plot of the weighted in- and out- degrees of C, 
Iq and O y c , respectively. The blue marker along the 
x-axis shows the periods in which the NIN03.4 (see 
e.g. |l4[) mean temperature exceeds 0.4°C, which is one 

of the criterions indicating an El-Niho event. The 
orange markers show the extended periods in which our 
analysis is influenced by the events because of the 365 
days integration of the cross covariance function. 



El-Niho events influence I and O in a very different 
way, having a rather peculiar phase locking. Strong El- 
Niho events weaken the links, yielding smaller micro- 
scopic contributions to both If, and O v c (white horizontal 
stripes in Figs. fUEkj) . just like Fig. [la] However, the in- 
links are significantly more vulnerable to El-Niho com- 
pared to the out-\mks. Moreover, before each major 
event of weakening /-field, there is a short, weaker re- 
sponse of weakening of the O-field. After the weakening 
epoch of both fields the O field recovers, and thereafter 
the / field recovers as well, leading to a full oscillation of 
both fields [15j . It is possible to summarize the situation 
by saying that the C nodes lose their autonomous power 
at the beginning of an El-Niho event, but thereafter, dur- 
ing the event, it recovers and becomes significantly more 
autonomous than the average. Despite the prominent 
difference during El-Niho between their dynamics, the / 
and O fields exhibit a remarkable static mirror similar- 
ity. Both fields are found to be highly asymmetrical with 
respect to the equator, at all times. We find that the in- 
links (going into C) mainly originate from the northern 
latitudes, while the out- links (out from C) mainly target 
towards the southern latitudes. 

In order to detect the weak response of O v c , we next 
consider specific "microscopic" contributions that have 
the lowest weight values, disregarding the rest of the field. 
This filters out noise that is not related to the observed 
slight weakening. We count (Fig. [3]) the number of nodes 
having microscopic contribution to the weighted in- and 
out- degrees under a weight threshold of 5, defined as I* 
and O* respectively. We will show (Fig. H]) that there 
is little sensitivity to this threshold value. Thus, a rise 
of O* (red) in Fig. [3] corresponds to nodes losing their 
dependence on C. Shortly after such events we see a 
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FIG. 3: Number of nodes outside C having sum of 
weights below 5 directing towards C (in), /*, or from C 
(out),0*. The blue and orange markers are described in 
the previous figure 
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FIG. 4: Cross covariance of the number of links below 
various thresholds going into and out from cluster C (in 
and out respectively, Fig. [3]). In the positive direction of 
the x-axis, the number of in- links lags in time after the 
out- links. 



significant rise of the I* (black) that corresponds to nodes 
losing their influence on C. As is clearly seen in Fig. [3j 
the /* and O* oscillations begin at the onset of El-Niho 
events. 

Next, we extract the time lag between the two coher- 
ent oscillations of I* and O* using the cross covariance 
function. In Fig. H] we show the cross covariance between 
the two curves of Fig. [3j for several threshold values. The 
peak of the cross covariance function shows that O lags 
after / in 100 ± 50 days [l6| . The sharp correlation be- 
tween the number of "weak" nodes in both directions (in 
and out) as well as the time lag is persistent up to a 
threshold value of 15. The other peaks of the cross corre- 
lation function are related to the quasi-periodicity of the 
El-Niho southern oscillations. 

Next we show in Figs. IHfcJbl the spatial distribution 
of If and Of during normal (non-El-Niho) times (We 
used the GMT tool [17[ to overlay the fields on maps). 
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FIG. 5: Spatial distribution of the weighted degrees of 
nodes (in and out) averaged over non-El-Niho times. 
The nodes of the network are shown by small dots. 





Our results also suggest the existence of a robust de- 
layed relationship between the inward and outward cou- 
5oo pii n g of ENB with the rest of the world. The emergence 
™ of the autonomous behavior is consistently following, af- 
™ ter typically 3 months, a short and weak episode of a de- 
2° q creased outward coupling. The two fields (outward and 
inward coupling) oscillate with a phase shift. 

We also find that ENB is forced more by the northern 
hemisphere than by the southern, and forces the south- 
ern hemisphere more than it forces the northern. Since 
the annual cycle in the two hemispheres is opposite, this 
north-south asymmetry might be related to the known 
(not yet fully understood) partial phase locking of the 
ENSO cycle with the annual cycle (see e.g. 2p|. Elj). 
Acknowledgements We thank the EU EPIWORK 
6 2 project, the Israel Science Foundation, DTRA, ONR, 
2° Q and Research project for a sustainable development of 
™ economic and social structure dependent on the environ- 
ooo merit of the eastern coast of Asia. 
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FIG. 6: Spatial distribution of the weighted degrees of 
nodes (in and out) averaged over El-Nino times. The 
nodes of the network are shown by small dots. The 
nodes belonging to C are shown on the weighted 
m-degree map as larger blue circles 



Our statement about the m-weighted degrees having a 
broader distribution by 10 to 20 percent at all times, gets 
its direct spatial interpretation. As seen in Fig. [5j the 
out- weighted degree spatial field is much more homoge- 
neous compared to the in- weighted degree which has a 
sharp minimum around ENB. This asymmetry between 
Jq and O v c becomes significantly more pronounced dur- 
ing El-Niho periods (see Fig. [6]). The cluster C is marked 
explicitly with large blue circles. In this sense, ENB gets 
significantly more autonomous during El-Niho. A full 
time dependent animation of these maps is available at 
[l8| . Note that the fields in Figs. [5] and [6] have been 
calculated without any correlation threshold. 

In summary, we have found a new dynamical pattern 
that reflects the coupling between the El-Niho basin, and 
the rest of the world. ENB becomes significantly more 
autonomous during El-Niho, losing a large fraction of its 
m-links, while still having o^-links. This kind of topol- 
ogy is reminiscent of pacemakers in network models [19| . 
The major impact of events inside ENB on world cli- 
mate on one hand, and the weakened correlations during 
El-Niho episodes on the other hand [3|-[5|, are thus not 
contradicting. In fact, the uni-directional interaction of 
ENB with large parts of the climate network might sug- 
gest the origin for its significant dynamical role in the 
global climate. 
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